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Abstract 



We introduce a new theoretical approach to dissipative quantum systems. By means of a contin- 
uous sequence of infinitesimal unitary transformations, we decouple the small quantum system 
that one is interested in from its thermodynamically large environment. This yields a trivial final 
transformed Hamiltonian. Dissipation enters through the observation that generically observ- 
ables "decay" completely under these unitary transformations, i.e. are completely transformed 
into other terms. As a nontrivial example the spin-boson model is discussed in some detail. For 
the super-Ohmic bath we obtain a very satisfactory description of short, intermediate and long 
time scales at small temperatures. This can be tested from the generalized Shiba-relation that 
is fulfilled within numerical errors. 
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1 Introduction 



1.1 Methods for dissipative quantum systems 

Dissipative quantum systems, i.e. quantum systems coupled to an environment, occur in many 
areas of physics and chemistry and have been studied extensively in the past using various 
methods. In their pioneering work Caldeira and Leggett Jl], |2| proposed a standard model of 
dissipative quantum systems describing a particle coupled to a thermodynamically large envi- 
ronment modelled by harmonic oscillators. The typical Hamiltonian for such a problem is of the 
form 

H = H S + H B + H SB . (1.1) 

H$ is the Hamiltonian of the system, Hb describes the environment and H$b is the coupling 
between system and environment. A recent review of methods in the field of dissipative quantum 
systems can be found in Ref. ||. 

Most approaches start off by discussing the time evolution of the reduced density matrix 
of the small quantum system. This is achieved by eliminating the bath degrees of freedom, 
for example by taking a partial trace or by integrating them out in a path integral framework 
using the Feynman- Vernon influence functional method [0]. This yields some effective action 
with interactions that are non-local in imaginary time ||. For a general dissipative quantum 
system some approximations become necessary in this effective action, for example the NIBA 
(Non-Interacting Blip Approximation) in the spin-boson model [|[. 

Another important approach is the description as a quantum Markov process. However, from 
a fundamental point of view this description should be derived from an underlying system plus 
bath Hamiltonian. It is well-known that this is only possible in certain limiting cases, for a 
careful discussion see e.g. ||. 

Our approach introduced in this work avoids this problem altogether by staying in a Hamil- 
tonian framework. A well-developed method for Hamiltonians is provided by renormalization 
theory. Three decades ago Wilson 0] studied the problem of a fixed source Hamiltonian using 
renormalization. This model is a simple dissipative quantum system. A nucleon with infinite 
mass, the fixed source, is coupled to it mesons. The nucleon contains only an internal degree of 
freedom, it can be a proton or a neutron. It is described by a two-level system. The two-level 
system is coupled linearly to two bosonic baths describing the two tt mesons. Wilson analysed 
this model using non-perturbative renormalization. A similar renormalization scheme has been 
applied to the problem of dissipative tunneling by Bray and Moore || and by Chakravarty ||. 
In the simplest case of dissipative tunneling one considers a particle in a double-well potential 
coupled to the bosonic excitations of the environment. The particle in the double-well potential 
can be approximated by a two-level system. In the review of Leggett et al. || the adiabatic 
renormalization scheme applied to the dissipative two-level system is also discussed. 

In a renormalization scheme the Hamiltonian of system plus bath is mapped to an effective 
Hamiltonian that contains no high energy modes but has the same low energy spectrum as the 
original Hamiltonian. The main problem is that the renormalized Hamiltonian still contains a 
large part of the bosonic modes of the environment, namely the slow and the resonant modes. 
Therefore it is rather difficult to calculate time-dependent correlation functions using this renor- 
malized Hamiltonian. In fact no example is known to the authors where such a renormalized 
Hamiltonian has been explicitly useful beyond the identification of the low-energy scale. Alter- 
natively one can use numerical renormalization |jlC| ]. But for dissipative quantum systems this 
method has (at least until today) a restricted range of applicability: One can only treat fermionic 
baths with this method. Therefore one has to find a mapping from the problem of interest with 



2 



a bosonic bath to a model with a fermionic bath. This is only possible if the bosonic bath has 
special spectral properties. Only recently this approach was used to calculate the dynamical 



properties of the dissipative two-level system [11], which can be mapped to the Kondo problem 
if one starts off with an Ohmic bath. 

Our approach used in this paper can be seen as an extension of traditional renormalization 
methods. By decoupling states with respect to their specific energy differences (and not treating 
high-energy modes only starting from the UV— cutoff), we are finally left with a simple effective 
Hamiltonian. This Hamiltonian contains only nearly resonant modes and can at zero temperature 
be solved without problems. Let us mention that a similar renormalization scheme has recently 
also been introduced by Glazek and Wilson in high-energy physics |T^ ]. Obviously our approach 
is fundamentally different from the usual methodology as we do not employ a reduced density 
matrix formalism. We rather obtain a final Hamiltonian where the small quantum system and 
the environment are decoupled. 

1.2 Flow equations for Hamiltonians 

Let us explain the approach used in this paper in more detail. Technically it is based on a 
continuous sequence of infinitesimal unitary transformations that is applied to the Hamiltonian. 
This technique has been proposed by Wegner to bring a given many-particle Hamiltonian in 
diagonal or block-diagonal form fl3|| . The continuous sequence of unitary transformations U(£) 
is labelled by a flow parameter £ with dimension (Energy) -2 . It corresponds to the energy 
difference that is just being decoupled, that is for small £ large energy differences are decoupled 
first before smaller energy differences are dealt with later for larger values of I. 

Applying such a transformation to a given Hamiltonian, this Hamiltonian (i.e. its parameters) 
become functions of £. Usually it is more convenient to work with a differential formulation 

— = [ V (£),H(£)}, H(£ = 0) = H (1.2) 
with an anti-hermitian generator rjil) related to the unitary transformation by 

V(i) = ^W)- (1-3) 

With a suitable choice of r](£) and additional approximations which will be described later, one 
is able to solve the differential equation ( |1.2[ ). Thereby one obtains a final Hamiltonian H(£ = 
oo). The necessary approximations are chosen so that the low energy spectrum of the initial 



Hamiltonian is well-reproduced by the final Hamiltonian H(£ = oo). In previous work [14, [y| 
we have already used this method to obtain the long-time behaviour of equilibrium correlation 
functions in the spin-boson model. In the present paper we show how the method can be 
improved to obtain information about equilibrium correlation functions also for intermediate 
(t A r ~ 1) and short times (t A r < 1). 

The generator rj{£) is chosen so that the initial Hamiltonian ( |1.1|) is mapped onto a final 
Hamiltonian which has the simple structure 

H{£ = oo) = Hsoo + H B . (1.4) 

The coupling Hsb is eliminated. Hsoo is a renormalized or effective Hamiltonian of the quantum 
system and Hb is the Hamiltonian of the environment. The latter remains unchanged since we 
assume that the bath is large compared to the system and therefore does not change its properties 
when it is coupled to the small system. Having such a procedure in mind two questions are 
immediate: 
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• At first sight it is quite unclear how a Hamiltonian like ( |1.4[) can describe dissipation. Since 
the small system is decoupled from the environment, exchange of energy between system 
and environment is no longer possible. To resolve this question, one has to understand 
how time-dependent correlation functions are calculated within our formalism. Usually 
dissipation manifests itself in a decay of correlation functions like {q(t)q(0)), where q is the 
coordinate of the particle coupled to the bath. We will later see that in our approach dissi- 
pation enters through the following simple fact: If one performs a unitary transformation, 
all the observables have to be transformed as well. In the case of a single particle, one 
has for example to transform the coordinate q. This means that one has to solve the flow 
equation 

g = M (1-5) 

Transforming an observable of the system, it will finally contain many environment oper- 
ators. Typically the transformed coordinate q(£ = oo) will be of the form q{£ = oo) = 
^2k c kQk, where q^ are coordinates of the bath oscillators and are some numbers. In 
order to calculate correlation functions we have to use the transformed observables and the 
transformed Hamiltonian. In this way it will be shown that the information on dissipation 
is contained in the unitary transformation that has been used to decouple the system from 
the environment. This will become clear in the special examples we discuss in Sects. 3 
and 4. 

• Which choice of r\ is appropriate for the problem we have in mind? The main idea is that 
rj is chosen in such a way that couplings between system and bath which correspond to 
large energy differences are eliminated first, couplings which correspond to small energy 
differences are eliminated later. This yields the fundamental separation of energy scales 
underlying renormalization theory. A suitable choice of rj was suggested by Wegner 
7](£) = [Hs(£) + Hb(£), Hsb(£)]- In the present paper we will use a different ansatz, which 
has some advantages to be discussed later. Nevertheless the couplings between system and 
bath corresponding to large energy differences are eliminated first. 

Another question that arises is whether one really needs infinitesimal unitary transformations. 
Now in principle it is possible to solve (|l.3|) and to obtain U(£ = oo) as an ^-ordered product 
of factors exp(rj(£)d£). Since the generator r](£) does not commute with itself for different values 
of £, it is in general not possible to evaluate this infinite product. Of course one could guess some 
single unitary transformation U right away. But U{£ = oo) constructed via the infinitesimal 
transformations r/(£) has the important energy scale separation automatically built in and this 
is difficult to ensure if one starts off with one single unitary transformation instead. In fact the 
unitary transformation constructed in the flow equation approach naturally contains "counter- 
terms" that eliminate singular induced interaction due to vanishing energy denominators. For 
the Anderson impurity model this has been discussed in Ref. 16]. 
Some advantages of this approach are: 

+ We work in a Hamiltonian framework. 

+ The procedure is non-perturbative as it relies on unitary transformations and has energy 
scale separation built in. Therefore one finds the correct low-energy scale in the effective 
Hamiltonian for £ — > oo. 
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+ The effective Hamiltonian for t — > oo is simple in the sense that only resonant couplings 
remain. This is the major advantage as compared to adiabatic renormalization or "poor 
man's scaling" . At low temperature our effective Hamiltonian becomes equivalent to a dis- 
sipative harmonic oscillator and equilibrium correlation functions are obtained easily. The 
main difference as compared to adiabatic renormalization is that the flow of the observables 
is not negligible and even very important for i corresponding to the low-energy scale. 

Now since one is usually interested in correlation functions of the system, it is sufficient to 
transform the observables of the system like for example the coordinate and the momentum of 
a particle. One problem is that apart from special cases like the dissipative harmonic oscillator, 
the flow equations for the observables are not closed and additional approximations become 
necessary. In our previous work [[lj], we have only discussed the long-time behaviour of 
equilibrium correlation functions of the dissipative two-level system. One goal of the present 
paper is to show how our method can be improved so that the whole frequency range can be 
treated. But the transformation of the observables causes also extra problems: 

- Since the flow equations for the observables are not closed, a special ansatz for the trans- 
formed observables is necessary. In this paper we use an ansatz that is linear in the bosonic 
modes. Such an ansatz works only for low temperatures, that means T has to be small 
compared to the typical low-energy scale of H$- At the moment we are not able to obtain 
reliable results for larger values of T without using traditional approximations of the kind 
discussed in Sect. 1.1 

- Sometimes one is interested in the dynamical behaviour of an observable in a situation, 
where the system has been prepared in a given initial state. Caldeira and Leggett § called 
this situation quantum tunnelling and distinguished it from quantum coherence, which is 
studied in this paper. In our approach we would have to transform the initial state as well. 
We leave this topic open for future work. 



1.3 Outline of the paper 

The plan of our paper is as follows. In Sect. 2 we develop the general method. A continuous 
unitary transformation is constructed together with a suitable approximation that decouples the 
system from the environment. We derive the flow equation for the Hamiltonian describing the 
quantum system and the flow equation for the spectral function of the bath. We compare our 
approach with adiabatic renormalization. As a main result we obtain a general expression for the 
behaviour of dynamical correlation functions at low frequencies. This result is a generalization of 
our former result for the dissipative two-level system We show that a dissipative behaviour 
of correlation functions is obtained only if the corresponding operators decay completely. The 
formalism presented in this section is very general. It can be applied to systems with arbitrary 
spectral functions. The only restriction is that the quantum system is coupled linearly to the 
environment. 

As our formalism is quite different from other approaches to dissipative quantum systems, 
we discuss the exact solution for the dissipative harmonic oscillator in detail in Sect. 3. Since 
the Hamiltonian is quadratic, this problem can be solved for any spectral function of the bath. 
We show explicitly how the correlation functions are calculated. Although we only reproduce 
well-known results [17], the reader is encouraged to study this section. It provides the basic 



ideas for understanding how the information on the dissipative behaviour of the system can be 
contained in the continuous unitary transformation. Furthermore some of the results are needed 
later on. 
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In Sect. 4 we discuss the dissipative two-level system. We calculate the Fourier transform 
of the equilibrium correlation function. We are able to treat a super-Ohmic environment with 
arbitrary spectral properties or an Ohmic bath with small coupling. We obtain differential 
equations that can be solved numerically or discussed analytically for the qualitative behaviour. 
The numerical results are presented in Sect. 5. We compare our results with the non-interacting 
blip approximation. Both methods yield similar results for frequencies near the renormalized 
tunneling frequency. In addition we find the correct universal algebraic long-time behaviour of 
equilibrium correlation functions at zero temperature as put forward by Sassetti and Weiss 
Their generalized Shiba-relation connecting intermediate and long time scales far beyond 
this "simple" universal long-time behaviour is also found to hold within numerical errors for the 
super-Ohmic bath. 

Finally Sect. 6 contains some conclusions and an outline of possible future developments. 



2 Flow equations: The general framework 



In this section we develop a general method to treat model Hamiltonians for dissipative quantum 
systems. The method is an improvement and extension of the one used in our previous work on 



the spin-boson problem [14, 15 1. The notation and the approximation we use are similar to our 
treatment of the Anderson single impurity model 16]. 



2.1 The model Hamiltonian 

We study the properties of a small quantum system coupled to environmental degrees of freedom. 
As usual the environment is modeled by a set of non-interacting bosonic degrees of freedom, i.e. 
harmonic oscillators. The small quantum system is described by a general Hamiltonian H$- It 
can for instance be the Hamiltonian of a particle in a confining potential or of a spin degree 
of freedom. In this section the quantum system is not specified. We assume that the small 
quantum system is coupled linearly to the bosonic degrees of freedom. The Hamiltonian of 
quantum system plus environment is given by 

H = Hs + A^Xkih + bD + ^Uk-.blh: . (2.1) 

k k 

The operators and b^ are usual bosonic creation and annihilation operators. In the sequel 
we will perform the limit where the number of bosonic degrees of freedom goes to infinity. The 
couplings Afc enter via the combination 

J(uj)=J2^^-^k)- (2-2) 
k 

J{uj) is assumed to be a continuous function of uj. A typical example is the so called Ohmic 
bath, where J(uj) = 2aujg c (u>/uj c ). g c {x) is a suitable cutoff function with <? c (0) = 1. Typical 
examples are g c (x) = exp(— x) or g c {x) = 0(1 — x). The final results should not depend on the 
special choice of the cutoff function but be universal functions of some low-energy scale. 

We treat the case of a general spectral function J(w), but we will need some additional 
conditions on J(uj) to be introduced below. The operator A has to be specified when one treats 
a special example. Concerning Hs we assume in a large part of the paper that it has a pure point 
spectrum. Generalizations to the case where the spectrum of Hs has continuous parts is possible 
if we assume that the ground state of Hs is separated by a finite gap from the continuous part 
of the spectrum. Although we are able to treat the case of a free particle coupled to a bosonic 
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bath exactly using our method, we exclude in the present paper the case where H$ has a purely 
continuous spectrum. We do not discuss e.g. the case of a particle moving in a periodic potential. 

The form of the Hamiltonian ( |2.1| ) contains the kind of problems discussed by Caldeira 
and Leggett There a particle in a confining potential and coupled to an environment was 
discussed. Caldeira and Leggett argued that as long as the moving particle perturbs each degree 
of freedom of the environment only weakly and the temperature is small, one can without loss of 
generality represent the environment by a set of harmonic oscillators. Furthermore they showed 
that apart from some pathological cases it is possible to choose the coupling between system and 
bath to be linear in the oscillator coordinates. Let us mention that recently Castro Neto and 



Caldeira [2C] proposed other models for dissipative quantum systems that do not fall into the 



class of models described by (|2.1| ). 

2.2 Flow equations for the Hamiltonian 

We now want to apply a continuous unitary transformation to the Hamiltonian Q2.1| ). Such 
a transformation is defined by a generator rj that depends on a continuous variable £. The 
continuous unitary transformation is defined by ( |1.2| ) with the initial condition H(0) = H in 
(|2.l| ). The Hamiltonian and consequently the coupling constants become functions of I. For the 
generator of the continuous unitary transformation we make the following ansatz 

r, = iJ2A k (b k + bl)+J2 B k( b k-bl)+J2 r lk, q -(bk + bl)(b q -bl): . (2.3) 

k k k,q 



A k and B k are elements of the algebra spanned by Hs(0) and A. In comparison to [14, 15], rj 
contains an additional term. This term is introduced to eliminate effective couplings between 
different bosonic modes. With this choice of rj one can construct the exact solution of the problem 
in the case of a harmonic oscillator or a free particle coupled to a bosonic bath, see Sect. 3. The 
commutator of rj with H is given by 



[ri,H] = i^u k A k {b k -b\) + ^u k B k {b k + b\) 

k k 

+i£)[A fc , H s ](b k + b{) + J2[Bk, H s ](b k - b{) 

k k 

+%^\ q {[A k ,A\) : {b k + b\){b q + W q ) : +iJ2*k[A k ,A](2n k + l) 

k,q k 

+i^2x q ([A k ,A] - ([A k ,A])) : ( bk + bl)(b q + bl) : 

k,q 

+ ]T \ q [B k , A] : (b k - b\){b q + b\):+Y J h[B k , A]+ 

k,q k 

+ ^2m, q Uk ■ (b k - b{)(b q - b\) : +J2vk,g^q ■ (h + b\){b q + b\) : 

k,q k,q 

+2j2vk,qXq(b k + bl)A . (2.4) 

k,q 

The expectation value ( . ) is defined by ( . ) = ^ ' cxp(- pHs) ' ^ or zero temperature it is simply 
the ground state expectation value. We now choose A k and B k such that the coupling of the 
system to the bath has always the form given in (|2.l| ) with an operator A that does not depend 
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on £. The ^-dependence is contained in the couplings X k . This yields the following conditions 
for A k and B k 

A k = —[B k ,H s ], (2.5) 

w k 

which eliminates couplings of system and environment in [rj,H] containing b k — b k , and 

u k B k + i[A k ,H s ] = -X k Af(u k ,£), (2.6) 

which guarantees that A does not depend on I. Furthermore we choose r\ k ^ q so that couplings 
between different bath modes in (|2.4j ) do not occur. To achieve this we need the following 
conditions 

Vk,qUk + Vq,kU q = 0, (2.7) 

which eliminates couplings of the form (b k — b\){b q — b q ), and 

Vk, q Uq + Vq,kUk + iXq([A k ,A})+i\ k ([A q ,A]) = 0, (2.8) 

which eliminates couplings of the form (b k + b k )(b q + b q ). When these conditions are satisfied 
two terms in [r/, H] remain, namely the term in the fourth line and the first term in the fifth line 
of fl2.4|). They contain a coupling of the system to two bosonic modes. We neglect these terms. 
This approximation has to be justified a posteriori. The corresponding couplings are not present 
in the initial Hamiltonian. A rough estimate of these terms can be obtained if one integrates the 
coefficients over I. 

A comparison of the terms in [rj, H] with the terms in the Hamiltonian yields the flow equa- 
tions for the couplings X k and for the Hamiltonian H$ of the system 

= -WH, e ) + 2 E ^, A> ( 2 - 9 ) 

= i h [A k , A] (2n k + l) + ^X k [B k ,A} + . (2.10) 

M k k 

The conditions on A k , B k and on r\ k ^ q can be solved. We obtain 

I /-00 

B k = X k f(u k ,£) / dre ir( ^ +i0 + ) (e irHs Ae~ iTHs + e~ iTHs Ae irHs ) , (2.11) 

^ Jo 

-I /-OO 

A k = --X k f(uj k ,£) / dre iT{ulk+i0+) (e irHs Ae~ iTHs - e~ iTHs Ae irHs ) , (2.12) 
2 Jo 



and 



LO q 



Vk, g = i(([A k ,A])X q + ([A q ,A})X k )- 2 

^k ~ ^q 

i X k X„Lo, 



2UJ 1~ -■■! 



f/H»*) I'" dre iT ^ +i0 +H[e irHs Ae- lTHs -e- irHs Ae lTHs ,A]) 

k~ UJ q V ^0 

+f{u q ,£) dTe iT ^ +i0 +^([e iTHs Ae- iTHs - e~ irHs Ae i/rHs , A])\ . (2.13) 
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With these choices of Ak and Bk, the neglected terms in (|2.4| ) are 

: {b k + bl){b q + b\) : / dr^ +i0 +\[^ Hs Ae- iTHs - e~ irHs Ae irHs , A] 

k,q ' J ° 

- ([e irHs Ae~ iTHs - e~ irHs Ae irHs , A]}) (2.14) 

+ o E W(^ k ,£) : (bk ~ b\){b q + b\): / dr^ +i0 ^ [e lrHs Ae~ iTHs + e-^ie^.i] 
1 , Jo 

Since ([e 4r " Hs ^4e _ *' r " Hs + e _ * r ^ s Ae*' r ^ s , A]) = 0, the expectation value of these terms vanishes. 
Neglecting these terms yields a good approximation if the integral over £ of these terms is small. 
This is the case if are small or if fall off rapidly as functions of £. We will come back to 
this point at the end of this subsection. 

The flow equations for \k and H$ can be rewritten using the ^-dependent spectral function 
of the bath, 

JH = ^A t (^-^). (2.15) 

k 

Notice J(u,0) = J(u>). Then the flow equations are 

— — = -2J(u,£)f(u,£) -2iJ(u,£) / dw uj2 _ (jJ i2 

(f(u>,£) dTe iT ^ +i0+) ([e iTHs Ae- iTHs - e~ irHs Ae irHs , A\) 

+f{io\l) dTe iT ^ l+i0+) {[e iTHs Ae- iTHs - e~ iTHs Ae irHs , A])^j (2.16) 



dH s i_ 
~df ~ ~2 



oo 



dujJ(uj,£)f(u,£) / dTe iT ^ +i0+) 
Jo 

({2n(uj) + l)[e irHs Ae~ irHs - e~ iTHs Ae irHs , A] 

-[e irHs Ae~ irHs + e~ irHs Ae irHs , A]+) (2.17) 



The integrand in the second term in ( [2.16 ) contains a singularity oc (a; — u>') 1 . It occurs due 



to the corresponding singularity in t?^. The integral has to be interpreted as a principal value 
integral. The singularity in r/fc „ occurs because we tried to choose r/ so that Hg remains form 
invariant. This has the advantage that the number of flow equations is small. The price we have 
to pay is that the singularity can cause problems. We will come back to this point at the end of 
the next subsection. 

Let us for a moment regard the trivial case where the commutator [A, H$] is initially a 
number c. Then [A, Hs] = c for all £. Furthermore the neglected terms in [77, H] vanish and the 
flow equation becomes exact. Integrating we obtain 

Hsoo - H S (0) = -A 2 [ (2.18) 
J u 

In order to obtain a finite result we need the condition 
/ <hj J ±^- < 00. (2.19) 
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This means that 



J(uj, 0) oc lo s , s > for small u. (2.20) 

Furthermore J(u>, 0) should contain a high-energy cutoff uj c . These conditions on J(u>, 0) are 
usually taken for granted. They can e.g. be motivated from a usual adiabatic renormalization 
scheme, see Appendix 1. The case s = 1 is the Ohmic bath J(u, 0) = 2aujg c (cj/uj c ). s > 1 is 
called the super-Ohmic bath, s < 1 the sub-Ohmic bath. 

A second point concerning ( |2.18| ) is that the term on the right-hand side is negative semi- 
definite. We have to introduce a corresponding counterterm in Hs(0) since otherwise Hsoo can 
become unbounded from below. This is well-known from perturbative renormalization. A similar 
term occurs in the general situation as well. Writing e lrHs Ae~ trHs in ( |2.17| ) as a series of iterated 
commutators, 

e iTH SAe -irH s =A + iT[H ^ A] _ ±t 2 [H s , [H S , A]] + ... , (2.21) 

the anticommutator in the last term in fl2,17| ) yields a contribution — A 2 J duj to Hsoo- One 
has to add a counterterm A 2 j du J< - U) J ) ^ to Hs(0). This is well-known from adiabatic renormal- 
ization (compare Appendix 1). 

2.3 Conditions on f(u,£) 

The form of the function f(u,£) has been left open until now. It is clear that f(cj,£) has to be 
chosen such that J(uj,£) tends to zero in the limit I — > oo. In this subsection we analyse the flow 
equations ( 2.16| , |2.17| ). The goal of this analysis is first to find an appropriate choice of f(u,' 



and second to see which conditions we need to be sure that the neglected terms in ( [2.14D are 
really small. We will show that f(uj,£) has to fulfil the conditions 

f(uj,£) ~ to 2 for large u, (2.22) 

f(uj,£) ~ (uj — (e n — e m )) 2 for uj near some resonance energy (e n — e m ). (2.23) 

The reader who is not interested in the details of the derivation of these conditions is encouraged 
to skip the remainder of this subsection. 



Using ( 2.19 ) and the fact that J(w,0) has to contain a high-energy cutoff co c , we make for 



the moment an ansatz 

J(tu,£) = Klu s 9c (uj/Cu(£)), (2.24) 

where g c (x) is some cutoff function with g c (0) = 1. Cj{£) is an ^-dependent cutoff, o)(0) = u c . 
We assume that Co{£) decays monotonously with £. Such an ansatz is always possible for £ = 0. 
It contains most interesting physical situations. If we neglect the second term in ( 2.16| ), such an 
ansatz is possible for £ > using 



g'{oj/u) uj div 
2g{uj/Cj) uj 2 

If we now use only the first two terms in the expansion ( 2.21| ) we obtain 



dH s f , t ,uj s+1 . ,duj (A 2 2n(uj) + l„ TT t . . 

/ dujK g > {uJ /uj)— — + \ L [[H S ,A],A] . (2.26) 



d£ / Q 2 d£ \ u 2uj 2 
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As long as uj is large compared to typical excitations of H$ the solution of this equation is 



essentially equivalent to perturbative renormalization. Using ( 2.17 ) instead of ( |2,26| ) yields a 
result that is essentially equivalent to adiabatic renormalization. The main point is that as long 
as uj is large compared to typical excitations of Hs, couplings to bosonic modes with energies 
above uj are integrated out. 

The question is now whether or not ( |2.24 ) can be used as a suitable approximation to ( p. 16 ). 
What is the effect of the second term in ( |2.16 )? To study this question we use again the expansion 
(|2.21| ) to write 



f 00 1 1 

i / dTe^+^de^Ae-^ 3 - e~ irHs Ae irHs , A]) = {[[H s , A], A]) A, + 0(-==). (2.27) 
Jo ' U A u 6 

Therefore we choose f(oj,£) oc uj 2 for large uj. This is possible using a cutoff function g c (x) = 
exp(— x 2 ), which yields 



(2.28) 



We choose uj = (2£)~z for I > uj~ 2 , since then f(oj,£) = uj 2 does not depend on I for I > uj~ 2 . 
This is the first condition ( |2.22| ) on f(uj,£) mentioned above. The flow equation for J(u,£) can 
now be written in the form 



dJ(u,£) to 2 du 



duj 



di 



2^—J(uj,£)-2J(uj,£) — ([[H s ,A],A])Kuj 



s-2 



dxx^ 



x — uj 2 /uj 2 



(2.29) 



In our ansatz ( [2,24 ) K was assumed to be constant. Now, at least for oj 2 /oj 2 <C 1 this equation 
can be interpreted as an equation for the ^-dependence of the coupling constant K, 



(2.30) 



Since < and ([[fig, .A], A]) < 0, K decays monotonously. Integrating this equation yields 



K{£) = K{0)- L + 



2r(%)([[H s ,A],A]),_ 



5(0) 



s-l 



Cj{£) s - 1 ) 



-1 



for s ^ 1 



(2.31) 



K(£) = K(0) 



2 r(i)([[ J ff 5 ,A],A])ln||| 



for s = 1. 



(2.32) 



Up to now we always used the expansion ( |2.21| ) to analyse the flow equations. This is possible 
for sufficiently large values of uj. It is clear that the result (|2.31 , 2.32 ) together with the ansatz 



(2.24) applies only if oj lies not in the region of typical excitation energies of Hs- If this is the 
case the expansion ( [2.21 ) is no longer possible. In fact one should not interpret the ^-dependence 
of K as a coupling constant renormalization but as a correction to the ^-dependence of J(oj,£). 

To study the behaviour for smaller values of uj, let us consider the case where Hs has a pure 
point spectrum. If one calculates expressions [e trHs Ae~ lrHs ] in the basis where Hs is diagonal, 
the r-integrals in ( 2,16j ) and ( |2.17| ) can be evaluated. ( |2.16j ) can be written in the form 



dJ(u,£) 



(hi 



e_xp{-Pe^ f , uj'J{uj\£) 

,m a m,n 7} J{0J,£) / Oh) /2 



n,m 
oo 



Z S 



UJ* — UJ 



x ( f(uj,£) j dTe iT(uJ+i0+) (e iT{en - em) - e^ 6 "" 6 ™)) 
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+ f(uj',£) / d re ^(^'+iO+)( e ir( £n - £m ) _ e -ir(e n -e m )j 
JO / 

-2J{u, £)f(u, Wm,n eM ~^ n) J( U , £) J ^' J3{JA) 



tM + ] (2.33) 



w 2 -(e n -e m ) 2 w' 2 - (e n - e m ) ; 

Here e n are the eigenvalues of fls, = 2~^n ex P( — /^ e «) anc ^ a n,m are the matrix elements of 
A in the eigenbasis of Hg- This equation (and similarly the flow of H$) contains divergencies 
oc (uj— e n +e m )~ l , where e n —e m is an excitation energy of H$- In order to avoid these divergencies, 
we demand that f(u>,£) behaves like (u> — e n + e m ) 2 for uj e n — e m . This is the second condition 
( 2.23| ) mentioned above. In this way f(u,£) becomes a function of the eigenvalues of Hg- If one 



neglects the second term in (|2.16| ) a first consequence of this choice of f(oJ,£) is that J(ui,£) does 
not change too much for small values of £ if u ~ e n — e m whereas it decays if u lies not in the 
vicinity of a resonance. 

Now for values of u large compared to typical excitation energies of Hg, the ansatz ( |2.24| ) 
still holds. But for values of u small compared to typical excitation energies of H$ } the second 
term on the right-hand side of ( 2.33j ) is positive and the derivative of J(u>,£) becomes positive 



for not too large values of £. In terms of the couplings A& this means that A& increases for small 
values of uj^. Such a behaviour can cause problems since the neglected terms ( [2.14 ) in [77, H] 



contain a factor A^Aq. In general one can expect that the neglected terms are unimportant for 
indices k and q for which uJk or oo q are not too small. Whether or not they become important 
for small ujk an d oj q depends on the behaviour of J(u,0) for small u. A quantitative analysis 
in the case of the dissipative two-level system in Sect. 5 will show that no problems occur in 
the super-Ohmic case J(uj, 0) oc uj s , s > 1, and in the Ohmic case J(u>, 0) = 2auj for sufficiently 
small a. 

The deeper reason for this infrared problem is our choice of rj. Our goal was to construct 
rj in such a way that off-diagonal matrix elements in the Hamiltonian are eliminated fast if 
the difference of the corresponding diagonal matrix elements is large. Due to the denominator 
(ufc — LOq)^ 1 in rjk^ it happens that the derivative with respect to £ of some off-diagonal matrix 
elements becomes large even when the difference of the corresponding diagonal matrix elements 
is small. 

To conclude our analysis of the general flow equations for the couplings and for the effective 
Hamiltonian (2.16) and (2.17), we have shown that for small values of £ the result of these 



1 



equations is more or less equivalent to adiabatic renomalization. (21) 2 plays the role of the new 
high-energy cutoff. When (2£)~ 2 is of the order of typical excitation energies of Hs, the couplings 
to higher energies in the total Hamiltonian H(£) are already exponentially small. When £ becomes 
larger, (2£)~z can no longer be interpreted as a new high-energy cutoff: For such larger values 
of £ the flow equations yield an approximate diagonalization of the renormalized Hamiltonian. 

In the next two sections the general method is applied to two special models, the dissipative 
harmonic oscillator and the dissipative two-level system (spin-boson model). Our method is 
exact for the dissipative harmonic oscillator. In the spin-boson model we restrict ourselves to 
the super-Ohmic case with arbitrary coupling and to the Ohmic case with small coupling for the 
reasons mentioned above. 
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2.4 Asymptotic behaviour 

Let us now discuss the situation where i^Soo has a pure point spectrum. Let us first assume 
that all its eigenvalues are non-degenerate. We denote the eigenvalues by e noo and we assume 
fnoo < e m co for n < m. The first point we have to investigate is the asymptotic behaviour of the 
eigenvalues of Hg{£)- By e n we denote the eigenvalues of Hs(£). As a first step we choose f(uJ,£) 
such that 

f e 

J(u, i) oc J nm exp(-2C nm / (w - |e m - e n \fd£') (2.34) 

JO 

for sufficiently large values of I. Such a choice of f(uj,£) is always possible. It is consistent with 
our discussion in the previous subsection, cf. ( |2.22| , |2.23j ). Neglecting the second term in ( p. 16 ) 
we write ( |2.17| ) in the form 

dH s i r du) ^A rdre^ +)((2nM + l)[e^Me-^-e-^M e ^,^ 

-[e irHs Ae~ irHs + e~ irHs Ae irHs , A] + ) . (2.35) 



dt 



The matrix elements of A in the basis where H$oo is diagonal are denoted as a nm . The asymptotic 
behaviour of e n is now obtained from the equation 



m J JO 



de n i ^ /" j 8J(uj, 

~d£ ~ 2 



(2n(cj) + l)[e ir ( <En_€m ) — e ~""( e n-em) _ J e «-(en-<!m) _j_ e -ir(e„-e m )' 

^ /■ gJ(avQ (2n(a;) + l)(e w -e m )+a; 



The main contribution to the integral over co for large values of £ comes from values uj 
I fnoo — e m oo | • Therefore we obtain 

jp CX 2 ^ Qnm^mnJnmC'nrn I duJ 6Xp( 1C nm I (ui |6 ra £m|) ^ ) 

« £ J Jo 



m 

u — \e n — e m \ 



((2n(o;) + l)( )+w). (2.37) 



n c m | 



This equation is similar to the equation that determines the asymptotic behaviour for A in the 
spin-boson problem |l4|]. It can be analysed as in the paper by Lenz and Wegner [21]. We obtain 



£n = £noo + and J(u!,£) oc ^= for u> = \e noo — e moo \. It should be clear that the asymptotic 
behaviour can be different if Hsoo has some degeneracies. 

2.5 Flow equations for observables 

Our goal is to calculate dynamical correlation functions for dissipative quantum systems. As 
already mentioned in the introduction, the simple form of the final Hamiltonian H{£ = oo) 
makes it easy to calculate expectation values. But it is clear that a given observable has to be 
transformed using the continuous unitary transformation. The transformation of an operator O 
takes the form 

% = h O] (2.38) 
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In the generic case a solution of this equation will not be possible. We assume that initially 
0{£ = 0) is an observable of the quantum system and not of the bath. We use the following 
linear ansatz for 0{£) 

= O s + ^2 °k+(bk + b\) + i J2 °k-(bk ~ b\) + ... (2.39) 

k k 

We neglect higher normal-ordered terms in the expression for O. The ansatz is similar to our 



previous work on the spin-boson problem 14]. But the continuous unitary transformation is 



different. The flow equations for Os and O k ± are 

dOs 

d£ 



J2[A k , O k+ ](2n k + 1) - iJ2i B k, Ofc-](2n fc + 1) + ^[B k , °k+\+ + £)[4fc, O fc _]+(2.40) 



i 

k 



and 

dO k - 



i[A k ,O s ]+2J2rik, q O q+ , (2.41) 
q 

-i[B k ,O s }-2j2ri q ,kO q -. (2.42) 
<? 

As mentioned above we assume O k+ (0) = O k -(0) = 0. In contrast to Ref. Q| the operators 
O k ± are coupled for different k. It is clear that it depends strongly on the initial operator O(0) 
whether or not the ansatz ( |2.39| ) will give good results. We shall see in Sect. 3 that the ansatz 
becomes exact in the case of the dissipative harmonic oscillator if O is e.g. the coordinate or the 
momentum of the particle. It is not exact if O is e.g. a power of the coordinate. Using O given 
in ( |2.39| ) we can now calculate equilibrium correlation functions of the form 

C o(1)o(2) (t) = (e itH 0^e- UH 0^) = Tr (p e ^ m e~ im O^) (2.43) 

where p e q is the equilibrium statistical operator. are observables of the system. It is useful 
to calculate the correlation function at £ = oo since in this limit the system is decoupled from 
the bath. The time dependence and the equilibrium average can be calculated. The final result 
is 

Comowit) = (e^O^e-^O^Jsoc 

k 

(e- to **(n fc + 1) + e^n*) 
+ i^ tHs -O k X l 00 e-^ m oo + ^s-o^e^o^ 

k 

{e- iWk \n k + 1) - e iWkt n k ). (2.44) 
Let us consider this general expression for the correlation function at zero temperature. Let 

(i) (i) (i) (i) 

Snm be the matrix elements of Soo and 0)± the matrix elements of in the eigenbasis 

of Hsoo- Then we obtain 



fG (l) (2) (1) (2) - (1) (2) - (1) (2) v (245) 

V^k+On^k+nO ' ^k-On^k-nO ^ LKy k+On^ k-nO ^ ^ k-On^ k+nO> ■ \^-^°J 
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The Fourier transform C (i) (2) (u>) has the form 

C WO^ ( U ) = O^QrPsnO 5 ^ + e °°° ~ e ™«=) + Y( O k+0n O< k+n0 + O k-0n°h-nJ0 

n k n 

+^i+OnOfio + *°i-On°2no)<^ + 60oo - enoo - w*) (2.46) 
The first term in this expression yields an oscillating behaviour of C (i) (2)(t) for large t. In 
dissipative quantum systems this term is generically not present. Generically one has C^oii = 0, 

(2) 

Quo = if e noo — eooo hes in the support of J(u), 0). In the examples in Sects. 3 and 4 we shall 
see that the asymptotic behaviour of J(u>,£) discussed above guarantees 0^ n = 0, O^ = 0. 
If S Q n ytz 0, the derivative d in ( |2.41| ) would contain terms that vanish oc £~* in the limit 
I — > oo. Consequently Ok+ would diverge. This shows that for £ — > oo the observable O in ( [2.39 ) 
decays completely and one obtains Os(£ — > oo) = 0. For small values of this yields 

^OWOWM ~ O^O^u) (2.47) 
+ E(°SooOSoo + °£oo°l-oo + *'°2oo°I-oo + ^ 00 og 00 )^ - 

A: 

For small values of ojk the matrix elements of Ak are (Ak)nm = i^kj ' {^k-,£)a n m- The a nm do not 
depend on k, they obey the condition Y^ r ^nrH^ m — H^ r a rm = a nm . The corresponding matrix 
elements of Bk contain an additional factor uok and can be neglected. Therefore for small Wfc, 
the ^-dependence of Ok+oo is contained in a proportionality factor Afc(0). The matrix elements 
of Ofc-oo contain an additional factor tvk and can be neglected. Thus we obtain 

Cowom M " O^ O%5(u) oc J( W , 0) (2.48) 

for small values of u. Apart from possible oscillations given by the first terms in ( p.46[ ) the 
correlation function C (i) (2)(t) therefore shows generically an algebraic relaxation oc i~ s-1 for 
long times. We will see that the generic behaviour of the correlation functions in the non- 
degenerate case is similar to the behaviour of correlation functions for the dissipative harmonic 
oscillator. This generic behaviour is obtained for operators Os that do not commute with for 
small u>k- The result ( |2.48{ ) remains valid if we allow for some degeneracies in the spectrum of 
H$- The main condition is that the ground-state of Hs must be unique. The spectrum of Hg 
can even contain a continuous part at higher energies. The general result ( |2.48j ) applies to the 
case of a particle in an arbitrary potential with bound states. The case of a periodic potential is 
not included. 



3 The dissipative harmonic oscillator 

The dissipative harmonic oscillator can be used as a good testing ground for any approxima- 
tion scheme in dissipative quantum physics since it is exactly solvable. The reason is that the 
Hamiltonian is quadratic. A survey of the properties of the dissipative harmonic oscillator can be 
found in the book of Weiss ||. The exact solution of this problem for a general spectral function 
J(co,0) was given by Haake et al. [17]. In this section we will show that our method is exact in 



the case of the dissipative harmonic oscillator. We are able to calculate all correlation functions. 
But our main reason for treating this model is to illustrate how our method works and how one 
is able to obtain a dissipative behaviour with a final Hamiltonian H(£ = oo) that contains no 
coupling of the quantum system to the bath. Furthermore we will introduce methods to solve 
the flow equations that will be useful later in our treatment of the dissipative two-state system. 
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3.1 Flow equations for the Hamiltonian 

Our general procedure is exact in this case since the neglected terms in [77, H] in (|2.4|) vanish 

2 2 

identically. If one chooses 77 as above and Hg = + ^A 2 !? 2 , the term is not changed, there 
is no mass renormalization in our approach. One only obtains a potential renormalization, which 
can be written as a renormalization of A. Instead we write Hs in the form 

H = Atfb + ]T X k (b + W)(b k + b\)+Y J " k b\b k + E (3.1) 

k k 

and use a slightly different 77 to preserve the form of the Hamiltonian. 
V = ^k\b-b^(b k + bl) + ^2v { k \b + b^(b k -bl) 

k k 

+ VkAh + b{)(b q - b q ) + r] b (b 2 - W 2 ) (3.2) 

k.q 

The commutator of 77 with H is easily calculated. 

[77^] = Y,^k )A + ^k ) ^k + 2r ]b X k + 2j2mM(b + b')(b k + bl) 

k q 

+ Y J {r, i k ) u k + r l { k )A )(b-V){b k -bl) 

k 



+ £>r,f% + (6* + b\)(b q + 



fc,9 



+ J3»/fc >9 Wfc(6* - b\){b q - b\) 



k,q 



+2 %'h(b + 6 f )(6 + b r ) + 2r/ fe A(6 2 + V 2 ) (3.3) 

We choose 

4 1) ^ + ? ?f A = ( 3 - 4 ) 
to eliminate terms in [77, H] containing (b — tf)(b k — bV), 

Vk,qUk + Vq,kUq = (3.5) 

to eliminate terms containing (b k — b^)(b q — b q ), 

Vk, q Uq + rj^ k io k + 277^ A, + 277WA& = (3.6) 
to eliminate terms containing (b k + b k )(b q + b q ), and 

7frA + J^ 2) A fc = (3.7) 



fe 

(.2 Lt2l 



to eliminate terms containing (6 — 6 t ). With these choices, the Hamiltonian remains form 

dH 

dt 



invariant. Comparing the coefficients on the left- and right-hand sides of =j§- = [77, H] and using 



X k = 0(1/ y N), we obtain the flow equations 

§ = "E4 2 'A t , (3.8) 
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f! = 25>fA t + 2X>«»A t , (3.9) 

k k 

£-°<*>- < 3 ' 10 > 

^ = {'A + 4 2) + 2 J> fc)(? A, + 2r/ 6 A fc . (3.11) 

<? 

We choose 

r?^ = -X k Af(u k ,£), (3.12) 

r?[ 2) = X k uj k f(uj k ,£), (3.13) 

% >9 = = _ 2A f gA ^ (/> fc ,^) + /W)), (3.14) 
With 

J(a;^) = ^A^( W -w fc ) (3.16) 
the flow equations can be written in the form 

— = 4y <kjLjf(u,e)j(u,e), (3.i7) 

^ = 2 y dw(w - A)/>, £) J( W) £), (3.18) 
^jM = 2 (u; 2 -A 2 )/>,W^) 



-8AJ(u;,f) J da/a/ j^'^ (/>, I) + /V, £)) 
J(lo, I) dA 



d£ 



(3.19) 



In the preceeding section we discussed how f(u>,£) = — (lo 2 — A 2 )f(uj,£) has to be chosen. Since 
now the operator A = b + W induces only transitions between next-neighbour states of H$, the 
only relevant excitation energy is A. Following the arguments for the general case the suitable 
choice of f(oo,£) is f(oo,£) = (u — A) 2 . We will not need to specify this special form of f(u,£), 
but it is often helpful to have a special form of f(w,£) in mind. We introduce 

R(z,£) = A !duj^4- (3-20) 
J z-uj 2 

and calculate its derivative with respect to £. This yields 
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Obviously 

A? + 4R(A* o ,£)-A 2 oo = (3-22) 
is a solution if R(z, oo) = 0. This yields a quadratic equation for A(£) 



A{£? = A^ - 4A{£) I ckoco-^^. (3.23) 



A^ = A^ + 4A 



Aqq can be determined from this equation for i = 

(3.24, 

We require for the reason of stability 

A >4 [ —J(u,£). (3.25) 
J u 

This condition reflects the fact that the initial Hamiltonian has to contain a counter-term 
A 2 J &.J(u,£). For A = 4 / ^J(u,£) we obtain Aoo = 0. This is the case of a free parti- 
cle. We discuss this special case at the end of this section. As an example that can be solved 
explicitly we let 

2 

J(u,0)= ]" a 2 . (3.26) 
If J(uj, 0) has this simple Drude-like behaviour we obtain 

Al = A 2 - 2^A 07 a . (3.27) 
The final solution is 



A ^M^ + y<^!>f_ 2lAo7 3 a . {3 ,s) 

3.2 Transformation of operators and correlation functions 

Since rj is bilinear in the bosonic operators, the operators Ok± in ( p.39| ) are simply numbers if 
Os = b or Os = & . Therefore the neglected terms in ( [2.39D vanish identically. We make the 
ansatz 

b{£) = P(£)b + P(£W + E ak ^ hk + E (3.29) 

k k 

The flow equation for b is 

§ = [»?, 6]. (3-30) 
It leads directly to the flow equations for the parameters 

f = 2 ^ - £ + E ^ + E + E ^ ( 3 - 31 ) 



A- 



2a% + E % «* + E c a * - E % + E «*> ( 3 - 32 ) 
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= Vk ] P ~ V^P + + ~ Y ^kUq + Y ^,qa q + Y + Y VKqUq, (3.33) 

q q q i 

= V^P + Vk ] P + Vk ] P - Vk ] P + Y ^ ka 1 + Yl VKqfXq ~ Y Vl^q + Y ^^q- ( 3M ) 

q q q q 

Introducing 

r±= VA^J if3±(3) (3 " 35) 
\ ±i 



we obtain 

Since initially r + (0) = r_(0) = 1 and Sfc i+ (0) = S£ ; _(0) = 0, we obtain r + = r_ and = Sfc 5 _. 
In the following the subscripts ± are dropped. The quantity 

+Y s l = l ( 3 - 39 ) 



r 2 



is conserved. This conservation law reflects the fact that [&(£), 6 f (^)] = 1 for all I. In analogy to 
R(z, £) we introduce the functions 

Sl (z,£)=Y SkXkV Y k - ( 3 -40) 
k z ~ u l 

si 



k & 

Calculating the derivative with respect to we obtain 

aS2 j^ = -4(r + 2S 1 (zJ))Y SkXkV ^ k { {U!k,i \ ( 3 -42) 

5*M . ^_ 2(r + 2Sl( , f))E ^W) 



(A 2 + 4fl(z, *) - z) V Wg^ri^ . (3.43) 



This shows that 



(r + 2g 1 (z,<)) 2 

52(z '^-A2 + 4i?(z,^-, (3 - 44) 
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does not depend on i. Since Afc(oo) = one has R(z, oo) = 0, Si(z, oo) = and since Sfc(0) = 
one has 5i(z,0) = 0, S^C^O) = 0. Therefore we obtain 

S 2 (z, oo) - -J^- = -(A 2 , + 4R(z, 0) - z)- 1 . (3.45) 

This allows us to determine r(oo) 2 and Sfc(oo) 2 . An immediate consequence of ( ggp is that 
r(oo) = if Aqo lies in the support of J(u). This is the case we are interested in. We introduce 
the function 

K(u) = J2sk(oo) 2 5(oo 2 -LJ 2 k ), (3.46) 

k 

which can be determined from S 2 {z, oo) using 
K(w) = -%S 2 {uj 2 -i0 + ,oo) 

7T 



77 \ A 2 - W 2 + iQ+ + 2AqP + 2vriA J(u>, 0) 

2A J(w,0) 

A 2 - u, 2 + 2A P (J ^S^) ) 2 + 4vr 2 A 2 J 2 („, 0) ' 



(3.47) 



P(.) denotes the principal value of the integral. The first term in the denominator vanishes for 
co = Aqq. K(uj) shows a maximum for some u> < A^. The shift of the maximum to a value 
below the resonance energy is due to damping. The behaviour of K(oj) for small u> is determined 
by the behaviour of J(u>) for small u>. From ( p.39 ) we obtain the sum rule 

POO 

2 / cIujujK(uj) = 1 (3.48) 
Jo 

If J(u>, 0) shows a Drude-like behaviour ( |3.26| ) we obtain 

K{ U ) = 2a 7 2 W + . 2 ) 

(A 2 ( 7 2 + a; 2 ) - 2to 7 3 A - w 2 ( 7 2 + w 2 )) + 4vr 2 Q 2 A 2 7 4 a; 2 



J_ (Ag-Ag )(7 2 + Ag M7 2 + ^ 2 ) 

" vr 7 ((A 2 , - a; 2 ) 2 (A 2 + 7 2 + c 2 - A 2 ) 2 + cV 2 (A 2 - A 2 ,) 2 ^ 2 + A 2 ,) 2 ' 

Correlation functions can now be calculated. We define 



(3.50) 



Cbb (t) d ^ (e im be- iHt b) T = C 6t6 t = (e im Ve- im V) T , (3.51) 

C 6t6 (t) = C bbt H f (3.52) 

where {.}t denotes the thermal equilibrium average. In the following we assume that Aoo lies in 
the support of J {to). Evaluating these correlation function at i = oo we obtain 

C bb (t) = J2 a fc (oc)a fc (oo)(e-^<& fc &t ) T + e^*(6^ fc ) T ). (3.53) 
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afc(oo) and ak(oo) can be expressed using Sfc(oo). The sum over k can be written as an integral 
over oo. This yields 



C bb {t) 



doo 



ujK{uj) f A 



00 



oo 



2 COS Lot 

e"/ T - 1 



Similarly we obtain 

loK(co) 



C&t&(*) 



C bb t(t) 



dco- 



duo 



loK(lo) 



w A 



2 cos u;i 4i sin uot 
e"/ T - 1 + e"l T - 1 



+ " _ 2 ) e --" 
lo A 



Ar 



a; 



+ 



LO 

An 



2 cos 



Ai sin 



+ 



An LO 

-^ + — + 2 

oo A 



-iuit 



(3.54) 



(3.55) 



(3.56) 



For T = the long-time dynamics is determined by the low-frequency dependence of J(cj,0). 
If we have J(lo, 0) oc lo s for small oo, then K(lo) oc u/ for small oo as well and the correlation 
functions C bb (t), C bb ](t) and C b \ h {t) behave as t~ s ~ l for large t. The same behaviour is obtained 
for C qq {t) = {C b \ b {t) + Cfefot ft) + 2Cj,fc(t)), whereas the velocity correlation function C OT (t) = 
^(2C fefe (t) - C 6t ° 6 (t) - C wt (t)) decays as r s " 3 . 



It is instructive to see how the effect of dissipation enters in the formulae. The main point is 
that in the denominator of K(lo) in the last expression in (p7|) the resonance oo = Aoo is damped 



by the second term. This term is due to the imaginary part of the denominator of ^(w 2 — «0 + , oo) 
in ( |3.45| ), i.e. the imaginary part of R(oo 2 — i0 + , 0). The final Hamiltonian contains no information 
on the dissipative behaviour of the system. It is simply a sum of the Hamiltonian of the bath 
and the renormalized Hamiltonian of the system, which is a Hamiltonian of a harmonic oscillator 
with frequency A^. The information on the dissipative behaviour is completely contained in 
the continuous unitary transformation and consequently in the transformed operators 6(oo) and 
fr^oo). These operators can be expressed as sums over bath operators alone, the corresponding 
term 05(00) in the general expression ( |2.39 ) vanishes since r(oo) = 0. 

Let us close this subsection with a remark concerning the spectrum of the Hamiltonian. It 
is well-known that the pure point spectrum of the original Hamiltonian H(£ = 0) contains only 
the ground state energy. The rest of the spectrum is absolutely continuous |]22| , |23[| . In contrast, 
we obtain a final Hamiltonian H(£ = 00) that has a pure point spectrum, which contains more 
eigenvalues, not only the ground state energy. At a first glance this seems to be a contradiction. 
But it is clear that for any finite value of i the Hamiltonian H (I) has only one eigenstate, namely 
the ground state. The argument by Arai |22| is sufficiently general and can be applied. The 
problem is that the excited states of H (I = 00) cannot be transformed back to I = since these 
states are completely embedded into the continuum of bath states. Strictly speaking the limit 
£ — > 00 does not commute with the thermodynamic limit. This was already observed in the first 
paper on flow equation for Hamiltonians by Wegner If one wants to study the problem of 
the spectrum of H rigorously, one has to stop the flow equations at some very large but finite 
value of i and treat J(oo,£) as in [22|. The results for the correlation functions remain the same. 



3.3 The free particle 

The case of a free particle coupled to a dissipative environment has been studied e.g. by Aslangul 
et al. 1 24]. It is described by the general case if Hsoo is the Hamiltonian of a free particle, i.e. 

Hs<x 



The treatment of this case is similar to the harmonic oscillator. The behaviour can 
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be deduced from the calculations above by taking the limit Aoo — ► 0. For the special example of 
a Drude-like behaviour ( |3.50| ) we obtain 

K{ ~ U) = ^u; 4 ( 7 2 +°J 2 - A 2 ) 2 + i 2 A 4 7 2 ' ( 3 - 57 ) 

Since now K{u) oc lo^ 1 , the correlation functions Cbb etc. are not well-defined. We have to 
calculate the velocity correlation function 

C vv {t) = (v(t)v(0)) T . (3.58) 

It can be obtained for the harmonic oscillator as a suitable combination of the correlation func- 
tions for b and b\ namely 

C vv (t) = ^C bb (t) - C m {t) - C Mt (t)). (3-59) 
A second point is that from ( |3.45| ) one obtains 

r(oo)- 2 = 1 - 4A J du^^- (3.60) 

For s < 2 one has r(oo) = 0, whereas r(oo) ^ for s > 2. Therefore C vv (t) contains an additional 
contribution 

= r -^lk B T (3.61) 

This yields 

C vv (t) = ± [ du^KM { ^4 + e~*A + r -^lk B T. (3.62) 
2m J \ e u i x — 1 / m 

For the above example fl3.26| ) we obtained K(u) oc lo" _1 , which yields C w (t) oc t~ 2 for T = 0. 
Consequently one finds C gq (t) oc hit for i — > oo. 

( |3.62| ) can be used for any general J(cj, 0). For = one has Ao = 4 J dw J ^'°) an d the 
general expression for if (a;) yields 

KM = ^ ^| . (3.03) 

Thus if J(u>,0) oc u> s for small oj, one has K(uj) oc u; _s if s < 2 and K(u) oc a/ -4 if s > 2. 
The velocity correlation function C vv (t) behaves like t~ 3+s for s < 2 and like i 1_s for s > 2 and 
T = 0. For finite temperature one finds C vv {t) oc i~ 2+s for s < 2 and C vv {t) —>■ r ^ k B T for 
s > 2. The factor mr(oo) -2 is often called a renormalized mass. Strictly speaking, however, 
neither in our approach nor in adiabatic renormalization there is mass renormalization. 

4 The dissipative two— level system 

The Hamiltonian of the dissipative two-level system is 

H = -^-a x + ^a z J2h(bl + b k )+J2"kblb k + E . (4.1) 
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This is of the general form ( |2.1| ) with 

H S = --(7 X + E , A = ^a z . (4.2) 

For a general discussion of various physical applications of this Hamiltonian we refer to the 
review of Leggett et al. || and to the book of Weiss ||. Our main interest in this section is 
to calculate the time-dependent equilibrium auto-correlation function of a z for this model. In 



14, 15] we already studied this model using flow equations, but we were only able to calculate 
the correlation function for low frequencies. Using our improved transformation we can now 
discuss the correlation function in the entire frequency range. 

4.1 Transformation of the Hamiltonian 

Due to the simplicity of Hg we are able to perform the transformation introduced for the general 
case in Sect. 2 without any additional approximations. Due to the transformation A and Eq 
become functions of I. The integrals in the general expressions for A k and B k can be calculated. 
The result is 

1 _A 

J 

and 



A k = --\ k f(u> k ,£) 2 _ A2 a y ( 4 - 3 ) 



B k = -\\ k f(u k ,l) 2 U \ 2 a z . (4.4) 
2 ojf, - A z 



Using these expressions we can derive rj ks 
l\ k \ a Auo a , , (3 A ( f(u k ,£) . 

Vk,q 



I*£*4 tanh 8* (ip4 + Ip4) (4.5) 

The flow equations for the spectral function of the bosonic bath and for the effective Hamiltonian 
of the quantum system are 

+2Atanh— J(u^) J du -^—^ ^-^-^ + -^-^ j , (4.6) 



^f = U duJ{u,£)f{u,£) f(2n{u) + V^jA^v* ~ ) ■ ( '•") 



This yields the flow equations for A and Eq. 

^ = -a J du,j(u,MM") + ^^r^ ( 4 - 8 ) 

— = ] du»,J M —^. (4.9) 

The second term in the flow equation for J(ui,£) has been neglected in [14, 15]. The effect of 



this term has been discussed qualitatively in Sect. 2. Let us now investigate what happens in 
the special case of the dissipative two-level system. The flow equations for the dissipative two- 
level system depend on the temperature. The reason is that we neglect higher normal ordered 
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interactions. We restrict our discussion to T <C A r . Since A is the only excitation energy of H$ 
in the case of a two-level system, the suitable choice for f(u,£) following the general discussion 
in Sect. 2.3 is 

f(u>,£) = (u - A) 2 . (4.10) 

With this choice the equations for A and Eq are free of divergencies. The only divergency 
is given by the denominator (u — lo')" 1 in the integrand in the second term in ( |4.6| ). If to is 
small compared to A, the second term is positive. For small £ typical values of to' are of the 
order Therefore the second term is of the order A£~z J(uj,£), whereas the first term is of 
the order A 2 J(lo,£). Thus the derivative of J(lo,£) is positive for small u and not too large £. 
Usually J(uj,£) has a maximum for some u of the order of A. But due to the fact that J(lo,£) 
increases with £ for small u, J(uj,£) can be very asymmetric with respect to its maximum. As a 
consequence, the derivative of A in ( |4.8| ) can change its sign. On the other hand the structure of 
the neglected terms ( [2.14 ) is similar to the integrand in the general flow equation of Hg ( [2.17 ). 



When the derivative of A changes its sign, the additional terms in [77, H] neglected so far are no 
longer small. The numerical solution of the flow equations in Sect. 5 shows that our method is 
applicable in the super-Ohmic case J(co, 0) oc co s with s > 1 for all coupling strengths, and for 
the Ohmic case J(co, 0) = 2auj with a ;$ 0.25. In the following analysis we restrict ourselves to 
these cases. 

Let us assume T = for the moment. Before we calculate the correlation function C(t) 
of the dissipative two-level system, we analyse the asymptotic behaviour of A and J(cu,£) for 
large £. If we neglect the second term in ( |4.6|) , we obtain exactly the same equations as in 
fll4| . The asymptotic behaviour of this equations is given by A(^) — Aoo oc £~^, J(u>,£) oc 
£~2 exp(— 2{ui — Aoo) 2 ^) |l4|, |2lJ. The second term does not alter the asymptotic behaviour of A 
as long as the above conditions (s > 1 or a ;$ 0.25 for the Ohmic case s = 1) are satisfied. As a 
consequence, the approximation 

^wr-J^^zg (411) 



is a good approximation to the flow equations (4.8) and (4.6) and can be used for £ = as a 
self-consistency condition for A^. 

4.2 The dynamical correlation function 

The main difference to our former approach [14, pl|] becomes apparent in the flow equations for 



observables. For cr z {£) we make the same ansatz as before, 

a z {£) = h(£)a z + a x ^ X k(£)(h + b\) . (4.12) 

k 

With this ansatz the flow equations for h and \k are 

* = -A£W2„ t + l)i<^ (4.13) 

k k 

~W = MAfc ^ + ^ x «^4 tanh — I^Fa^ + ^ J • (414) 



The second term in ( 4.14| ) was not taken into account in [14, 15]. When this term is not present 



Xk{£) diverges if 00k = Aqq. The divergence is only logarithmic, but it is clearly unphysical. This 
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was the main reason why we were not able to obtain reasonable results for correlation functions 
in a frequency range around A^ in Refs. (IJ], |l5[ . Due to the second term this divergence is now 
smeared out. This is the lesson from the dissipative harmonic oscillator, compare Sect. 3. As a 
consequence the Fourier transform of the correlation function 

= /i 2 (cx>)cos(A 00 i) + ^XfcM(2n fe + 1) cos(w fc t) (4.15) 

k 

is well-defined. 

Again due to normal-ordering, the flow equations for a z depend on the temperature. The 
neglected terms in ( [4 . 1 2| ) are normal-ordered and do not contribute if we calculate the simple 
average (cr z ). But in the correlation function ( 4.15| ) we have to calculate an average of a product 



of two normal-ordered operators. Therefore the neglected terms in ( f4.12| ) may be important. 
This happens when the temperature becomes of the order of Aqo. 

Let us now calculate the correlation function for low temperature T <C A r . If one neglects 



the second term in ( |4.14|) , the quantity h + YlkX-k 1S conserved. Due to the second term in 



( f4.14| ) this sum rule is only asymptotically fulfilled. We obtain 
± (h 2 + V Xl) = 2A V ( IpA + (4.16) 

Using the asymptotic behaviour of Afc one shows that the right-hand side falls of like £~ 2 , which 
shows that h 2 + J2k xl ^ s asymptotically constant. As a consequence we obtain h — ► if Aoo lies 
in the support of J(u, 0). Using ( 4.13j ) one shows that h(£) falls of like t~ * . Thus the asymptotic 



behaviour of h is not affected by the second term on the right-hand side of ( 4.14 ). This term is 
only important for uj^ A and I « A^ 2 . The fact that h 2 + x\ is asymptotically constant 
can be explained if one compares the flow equations for h and Xk with the flow equations for r 
and Sfc, (3.37) and ( |3.38| ). The latter can be written explicitly in the form 



dr - -2Y,XkSk§^ 2 VA^k (4.17) 



ds k tt^kJ) pr— .v- AfcAgA^u^ f f(oJ k J) , f(w q ,i) \ fA , a , 

~t = ^^T^V^ + 4 L «« + ^ J • ( 4 - 18 ) 

Afc in the case of the dissipative two-level system corresponds to 2A& in the case of the dissipative 
harmonic oscillator due to the different definition of the couplings in the Hamiltonian. Further- 



more h corresponds to r and Xk corresponds to Sfe. The difference between (|4.13j) and (4.17) is 



that the factor \J Alu^ in ((4.17|) is replaced by A in ( 4.13; ). Similarly \J~KuJk in the first term in 



( 4.181 ") is replaced by A in fl4.14|) and ^JuTj^J^ in the second term is replaced by oj q . Since for large 



t the couplings Afc differ only significantly from zero if ujk lies near A, the difference between the 
two sets of equations becomes small for large I. Therefore to obtain an expression for the Fourier 
transform of the correlation function C(t), we proceed as in the case of the harmonic oscillator. 
We introduce the functions 

S 2 (z,£)=^-^, (4-19) 
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S 1 (z,l)=y V ^ Xk 2 Xk , (4.20) 

S (z,t)=J2^\- (4.21) 
The conserved quantity ( |3.44| ) in the case of the dissipative harmonic oscillator corresponds to 

^(M)- /t^f. (4.22) 
A^ — z + Aoo(z, «) 

for the dissipative two-level system. It is clear that this quantity is not conserved. But since for 
large £ the flow equations for the harmonic oscillator and for the dissipative two-level system are 
similar, this quantity will be approximately constant. Using the asymptotic behaviour of one 
shows that the derivative of (4.22) behaves like £~ 2 . The asymptotic behaviour of the relevant 
quantities, \k an d h is much slower. 

The physical reason for this approximate equivalence of the dissipative harmonic oscillator 



and the dissipative two-level system for large £ is the following simple observation: From ( |3.55 ) 
one obtains the mean occupation number for the harmonic oscillator at zero temperature 

+ (4 .23) 

Let us consider a spectral function J(uS) that is strongly peaked for co = Aq. Then K{uj) is 
strongly peaked as well and the integral ( [4.23 ) becomes small. Since for large £ the spectral 



function J(u>,£) is strongly peaked, the dissipative harmonic oscillator behaves for large £ like a 
two-level system: The higher states of the harmonic oscillator are not occupied. It is clear that 
this argument only holds for low temperatures T <C A r . 
Since (h + Si) — > for £ — ► oo, we obtain 

S 2 {z, oo) « S 2 {z, £ ) - 2 . r , . , . 4.24 

A z (£ ) - z + A(£ )S (z,£ ) 



for a sufficiently large value £q. We let £q = (2XA OQ )^ . The right-hand side of (4.24) can 
be obtained numerically by integrating the flow equations, which will be done in Sect. 5. The 
one-sided Fourier transform 

CM = ^ xl(oo) coth(^) 5{u - u k ) (4.25) 

k 

of the correlation function 

C{t) = I oIujC{oj) exp(iwt) (4.26) 
Jo 

can be obtained from S-z{z, oo) using C(uj) = —^-QS2(oj 2 — i0+, oo). Notice our unusual normal- 
ization condition as compared to the literature following from the initial condition C(t = 0) = 1 

oo 

duj C{co) = 1. (4.27) 



o 



The analytical solution of the set of differential equations leading to C(uS) is not possible. A 
qualitative impression can be obtained by solving the linearized differential equations for Afe 
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and Xk U P to some value £q = ^AA^) -2 with A = 0.2 . . . 0.5. For smaller values of A the 
nonlinearities are too important and for larger values the asymptotic conserved quantity ( |4.22| ) 
cannot be used. In order to solve the differential equations ( f4.6| ) and ( [4.14] ) one stills has to make 
the additional approximation to replace A{£) and h{£) on the right-hand side by (1 + A)Aoo 
and h = 1, resp. This is reasonably good by comparison with the numerical solution. Putting 
everything together one finds the following result for the correlation function 

o M . %W*(fi±%=^=S (-8) 

^exp(-2(^-(l + A)A oo ) 2 £ ) \ 

((1 + A) 2 A2 o - a, 2 + (1 + X)A oo RS (^ - *0+,£ )) 2 + + A) 2 A 2 . J{u,i Q y) 

Only the two most important terms from ( |4.24] ) are used here, the first coming from the solution 
of the linearized differential equations and the second term coming from the conserved quantity 
proportional to h(£o) 2 . Eq. ( 4.28| ) reflects the qualitative features of the numerical solution 



for uj — > 0, uj ~ A r and u 3> A r . Due to all the necessary approximations it is not a very good 
approximation to the numerical solution beyond such qualitative features that are independent 
of the arbitrary parameter A. We have therefore made no attempt to improve this analytical 
result by e.g. taking more terms into account. 

5 Numerical results 

In order to obtain quantitative results for the correlation function C(iv), we have to integrate 
the flow equations numerically. The strategy is obvious: First of all we numerically integrate 



the differential equations up to some sufficiently large value £$. This gives us S2(z,£q) in ( (4.24|) . 
Then we add the asymptotic conserved quantity to obtain 52(^,00) and finally C(lu). 

How this is done explicitly will be explained in the next subsection. In particular it has to 
be established that the final result is independent of £q. Subsection 5.2 contains quantitative 
results for the correlation function for Ohmic and super-Ohmic baths. All calculations are done 
for zero temperature but can easily be extended to T < A^. Such effects of nonzero but small 
temperature will be briefly explained in the text. The unrenormalized tunneling frequency Ao 
always defines the energy scale with Ao = 1. 

5.1 Technicalities and tests 

We solve the flow equations for the imaginary parts of the functions So(z, £), S\{z, £) and S^fo 
defined in the following manner 

Ji(u,£) = --5Si(uj 2 -i0+,£), i = 0,1,2. (5.1) 

7T 

For convenience we write Jo(u>, £) = f J(ui, £) in this section. From Ji(ui, £) we can also reconstruct 
the real part of Si(uj,£) 

^S l (u; 2 -iO + ,£) = P [ du/ ^^-W- - ( 5 - 2 ) 

J UJ Z — UJ /Z 

The complete set of differential equations that has to be solved numerically for the spin-boson 
problem is 

^ = "A fducoth^M.J)^ (5.3) 
at I 2 uj + A 
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d£ 



(hU ' xA) = -2(u- A) 2 J (uJ,£) + 2Atanh^J {ujJ) J duj'J {uj',£)I{uj,uj',£) (5.4) 
= - j dujVKJcoth^^ujJ)^-^ (5.5) 



^ = -(-A^fl+A^^fl^ + ^^fl (5.0, 



+Atanh^-(ji(cj,£) J du' J (lo',£)I(uj,uj',£) 

+Jq(u>,£) J duj' ^—J\(uj' ,£)I(u!,uj' 



where 



r, i „s J fu;-A(£) co'-A(£)\ 

= ri/2 + -r-irri; ■ ( 5 - 7 ) 



uj 2 - uj' 2 \uo + A{£) u' + A(^) 
J2(oj,£) can be obtained via the relation 

^,0=1^. (5.8, 

After integrating these equations up to £q, we add the approximate conserved quantity from 
( pi~2ll ) and obtain C(cj). 

Now the differential equations for Jq(uj,£) and J\{uo,£) have to be discretized for certain 
values uJi and can only then be solved numerically. Typically we have used about 200 bath 
modes. Obviously it makes no sense to use equidistant values of uji, rather one has to sample the 
vicinity of A^ much denser: For large £ the curves Jq(uj,£) and J\(uj,£) are only nonvanishing 
close to Aqq. This is reminiscent of numerical renormalization, however, there the vicinity of 
uj = is sampled denser although the main spectral weight might be somewhere else (compare 



e.g. Ref. [11]). In all cases we had at least a minimum of 10 bath modes within the half width 
of the final spectral function Jo(u;,^o)- 

The quality of the numerical routines and this discretization procedure can be tested conve- 
niently for the dissipative harmonic oscillator discussed in Sect. 3. The set of differential flow 
equations is very similar to (|5.3|)-( ^6[ ), also the conserved quantity has the same structure. We 
have compared the numerical result for the function K {uj) with the exact solution from ( |3.47 ) 



for various values of £$ = ^AAoo) . The agreement always turned out to be excellent. A 
typical result for a Drude-like spectral function ( 3.26| ) can e.g. be seen in Fig. 1. The maximum 



deviation of the curves is less than 2% for all values of uj. 

For the dissipative harmonic oscillator the conserved quantity is an exactly conserved quantity 
as we have seen in Sect. 3. This situation was different for the spin-boson problem where it was 
only an asymptotic conserved quantity. Hence a second important test for the consistency of 
our approach is to check that the final result C{uj) in the spin-boson model is independent 
of £q = (2AAoo) _2 if A is sufficiently small. 

We have e.g. investigated this for the Ohmic bath with the parameters chosen as in Figs. 2a- 
c. In every diagram the curve C(uj,£q) is the correlation function without using the conserved 
quantity and C{uj) is the final result after adding the approximate conserved quantity. One 
sees that for smaller A more and more spectral weight is contained in the numerical solution 
and less in the conserved quantity. However, the final results C{uo) for the different values of A 

agree very well as can be seen in Fig. 2d. The spectral function Jq{uj,£q) = f J(uj,£q) becomes 
more and more peaked around A^ for smaller A as can also be seen in Fig. 2a-c. Hence 
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the physical reasoning underlying the approximate conserved quantity (the harmonic oscillator 
being equivalent to a spin-boson model if Jq{uj,Iq) is strongly peaked around Aqo, compare the 
discussion following ( 4.23|) ) is confirmed very well. Similar tests have been made for the other 



sets of parameters discussed in subsection 5.2 as well always leading to the same conclusion. As 
a suitable compromise between computing time and the fact that we only have an asymptotic 
conserved quantity we have used A = 0.1 for all quantitative results in subsection 5.2. 

A final test for the consistency of our approximations is provided by the normalization condi- 
tion Jq 00 C{uj) dio = 1 derived from the initial condition C{t = 0) = 1. Remember that according 
to ( |4 . 1 6| ) this is not exactly fulfilled. However, for all curves in Fig. 2 and also in the next 
subsection the sum rule is fulfilled with an error of less than 1%. This agreement can be consid- 
ered excellent and also gives support to the ansatz for the transformed observable a z according 



to ( 4.12 ) as very little spectral weight is lost. 



5.2 Quantitative results 

Ohmic bath: 

Fig. 3 contains the correlation function C(oj) for an Ohmic bath J {to) = 2aujQ{uj c — oj) with 
various values of a. C{uj) is plotted as a function of uj/A r with the the low-energy scale A r . A r 
is here used as a fit parameter to identify the positions of the various peaks (therefore the areas 
under the curves are necessarily different). The half width of the curves gives the dimensionless 
Q-factor of the damped oscillations. For a = one has only a <5-peak corresponding to undamped 
oscillations and in Fig. 3 one observes how this peak decays as a function of increasing coupling a. 

Unfortunately, for larger values of a than contained in Fig. 3 we did no longer observe that 
the final result for C{u)) was independent of (.q. This problem can be traced back to the small 
energy denominators in the generator in (|4.5| ). We have already discussed this possible 
problem in Sect. 2.3 and we refer to this section for more information. Hence using the present 
approximations our approach is limited to small couplings a ;$ 0.25. 

Now let us analyse our results for such small couplings in more detail. In Fig. 4 we compare 
the flow equation result with the correlation function obtained by the NIBA-approximation (see 
e.g. Ref. ||). The quantitative agreement for intermediate time scales is very good. Differences 
are apparent in the long-time behaviour, i.e. the behaviour for small oj. This can be more 
clearly seen in the response function S{oj) = C '{uj) / {tt uj) in Fig. 4: The NIBA-result diverges for 
u> — > which is an indication of the wrong long-time behaviour C{t) oc t _2 ( 1_a ) ) whereas S{0) is 
finite in the flow equation framework. As expected the flow equations therefore yield the correct 
universal long-time behaviour C{t) oc t~ 2 . Hence we agree with the NIBA where the NIBA can 
probably be trusted and we disagree where the NIBA is known to be wrong. 

A very sensitive test is provided by the Shiba-relation [19] generalized to the spin-boson 
model in Ref. [18]. For an Ohmic bath this reads 



Hm £M 1 2a{2 Xo ) 2 . (5.9) 

The static susceptibility xo is extracted with a Kramer 's-Kronig relation and a fluctuation- 
dissipation theorem from the correlation function C{oj) 

tt -;r^ (5.10) 

2 ./n U 
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Notice that our normalization ([4.2 
which makes up for a factor n in ( 
and small a. 




du)C(u) 



1 differs from the definition in Ref. fll 



In Table 1 relation (5.9) is tested for various parameters 



a 




lim^o C(u)/u 


2a(2 Xo ) 2 


%error 


0.01 


40.0 


0.0209 


0.0212 


2% 


0.01 


100.0 


0.0213 


0.0222 


2% 


0.025 


40.0 


0.0644 


0.0625 


3% 


0.025 


100.0 


0.0682 


0.0660 


3% 


0.05 


40.0 


0.18 


0.16 


10% 


0.05 


100.0 


0.20 


0.18 


10% 


0.1 


40.0 


0.94 


0.69 


25% 


0.1 


100.0 


0.94 


0.69 


25% 



Table 1: Generalized Shiba-relation for an Ohmic bath as found from the numerical solution of 
the flow equations (Ao = 1). The typical numerical error of both numbers obtained numerically 
is estimated as 3%. 

One observes that for very small couplings the Shiba-relation is fulfilled within numerical 
errors. However, the deviations grow with increasing a which can again be traced back to the 



small energy denominator problem for the Ohmic bath in (4.5). The typical effect of this is 
that lim^^o C(u)/u comes out too large. We leave this problem for future investigation and 
note that for the super-Ohmic bath investigated below the generalized Shiba-relation turns out 
to be fulfilled with very good accuracy in the whole parameter range. Since the error in the 
Ohmic Shiba-relation can be traced back to small frequencies which carry very little spectral 
weight in the correlation function, this error does have little implications for the accuracy of the 
correlation function (for the NIBA the error is in fact infinite, still the NIBA probably makes 
sense on intermediate time scales for small couplings). 

The low-energy scale of the spin-boson model with Ohmic coupling is set by 

A r = cA (^Y~ a , (5.11) 



where c is some constant of order 1 as has already been shown in our previous paper [14]. 
This theoretical prediction was found to be well confirmed from the numerical results, where 
for small a the energy scale A r plays the role of a parameter determined from the maximum 
of C(u). One expects to find universal correlation functions as a function of uj/A r for u -C u> c . 
In Fig. 5 the response function S(u) is plotted as a function of this rescaled energy. Universality 
is confirmed with excellent accuracy. 

Let us mention that the curves in Fig. 5 can be compared directly with numerical renormal- 



ization results of Costi et al. [11]. Qualitatively the curves look similar for these parameters, 
however, the maximum in the NRG framework is considerably lower: S(A r )/S(0) ~ 4.0 from 
NRG whereas we find S(A r )/S(0) ~ 6.0 using flow equations. We believe this problem can be 
traced back to problems inherent to numerical renormalization approaches for small couplings a: 
One has to resolve a sharp peak in the correlation function where due to logarithmic discretiza- 
tion very few bath modes lie. From Fig. 4 it is apparent that discretization is no problem in the 
flow equation framework. 
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A final remark about small temperatures T < A^. The main effect is that the correla- 
tion functions C(u) in Fig. 3 have to be multiplied with a factor coth(/?o;/2), for a qualita- 
tive discussion the effect of nonzero temperature in Eqs. (|5.3j)-( |5"l)| ) is negligible. Therefore 
lim^-,.0 C(oS) = const. 7^ for nonzero temperature leading to the expected exponential long- 
time decay of C(t). 



Super-Ohmic bath: 

For a super-Ohmic bath the coupling function can be parametrized as 
J(w) = K x - s u s Q{u c - u), 8>1 (5.12) 
with a coupling constant K with dimension energy. The low-energy scale A r is set by 

A r = cA exp(-^ TJ (^)*- 1 ) (5.13) 

where c is a constant of order 1. For the numerical calculations we have used a super-Ohmic 
baths with s = 2, 3. Extension to other values of s > 1 is unproblematic and yields similar 
results. 

In Fig. 6 the correlation function C(u>) is plotted as calculated from flow equations or the 
NIBA. The qualitative agreement is reasonable for intermediate time scales. Anyway it is prob- 
lematic to justify the NIBA-approximation for super-Ohmic baths since the blips cannot be 
considered as a dilute gas ||. As in the Ohmic case, an important difference shows up in the 
long-time behaviour. Since lim^^o = const. 7^ in the NIBA (compare Fig. 6) one ob- 

tains an exponential long-time decay of C(t), whereas the flow equation result is C{oS) oc u> s for 
small u>, hence C{t) oc t _s_1 for long times. 

Again a sensitive test for the long-time behaviour is provided by the generalized Shiba- 
relation also put forward for super-Ohmic baths in Ref. |lq] . It reads 

Hm^ = K^(2 X0 f. (5.14) 

W— >0 UJ S 

In Table 2 this relation is tested for various parameters. 

One observes that the generalized Shiba-relation holds within numerical errors]]] Considering 
the fact that the Shiba-relation is a highly nontrivial universal property connecting intermediate 
and long time scales far beyond the "simple" universal C(t) oc behaviour, this result can 

be considered very satisfactory. This supports the conclusion that our ansatz for the flow of the 
Hamiltonian and the observables already contains all the "relevant" physics and is probably very 
close to the exact solution. Let us also remark that for super-Ohmic baths the Shiba-relation 
cannot be tested with numerical renormalization group methods. 

Correlation functions for various parameters of a super-Ohmic bath are plotted in Fig. 7. 
For increasing coupling corresponding to smaller values of K, the damping becomes stronger 
and stronger resulting in an increasing Q-factor. (Of course no phase transition occurs for some 
critical coupling as it does in the Ohmic case for a c = 1.) 

A final remark about nonzero but small temperature T < A r . Again the main modification is 
that the zero-temperature curves are multiplied with the factor coth(/3u;/2). However, for s > 1 

1 Notice that the deviations are systematic. But that cannot be interpreted with certainty due to our present 
numerical error of the used computer routines. 
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0.166 


0.172 


4% 


2 


20.0 
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2.19 
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5% 


3 


40.0 
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0.000710 
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20.0 
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0.00380 


0.00399 


5% 


3 


10.0 


20.0 
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0.0628 
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Table 2: Generalized Shiba-relation for a super-Ohmic bath as found from the numerical solution 
of the flow equations (Ao = 1). The typical numerical error of both numbers obtained numerically 
is estimated as 3%. 

this still leads to an algebraic long-time decay of C(t) proportional to t~ s . That is for very small 
temperature one still has an algebraic decay for super-Ohmic baths and not an exponential one. 
This might be surprising at first sight, however, for the exactly solvable dissipative harmonic 
oscillator one makes the same observation.^ 

6 Conclusions 

This has been a long and in parts probably tedious paper for the reader. However, the approach 
put forward here is quite complementary to common wisdom in dissipative quantum systems 
and required more explanations and justifications than usually necessary. 

Let us sum up the main results of this paper. We have investigated dissipative quantum 
systems with the specific examples of the dissipative harmonic oscillator and the spin-boson 
model using infinitesimal unitary transformations. The Hamiltonian for such problems is of the 
form 

H = H S + H B + H SB , (6.1) 

where H$ is the small system coupled via H$b to the environment described by the bath Hamil- 
tonian Hb- With a continuous sequence of unitary transformations, a unitarily equivalent Hamil- 
tonian has been found where system and bath are decoupled 

Hqo = Hsoo + Hb- (6.2) 

In order to do this some approximations that neglect higher ("irrelevant") normal-ordered in- 
teractions generated by the unitary transformations have been used for the spin-boson model. 
The quality of these approximations seems to be quite good as we have seen in this paper (see 
also below). Notice that our approach is systematic in the sense that higher order terms can 
successively be taken into account. 

2 It is not necessarily true that the long-time behaviour for nonzero temperature is an exponential decay set 
by the smallest Matsubara frequency as sometimes claimed in the literature. 



32 



What is perhaps most surprising at first sight is that such a program going from \6. \ ) to ( ffi.ij ) 
using unitary transformations can be carried out at all. The key question is where dissipation 
can enter in a description like (|T^). The answer lies hidden in the fact that the observables 
have to be transformed as well under the unitary transformation. If the discrete eigenstates 
of Hsoo are embedded in the continuum of bath states, we have found that generically these 
observables "decay" completely under the sequence of unitary transformations. That means the 
term describing the original observable with respect to ( |6.lD has vanished completely and been 
transformed into other terms with respect to the transformed Hamiltonian ( |6.2| ). These new 
terms show typical dissipative behaviour (e.g. decay of correlations in time) and in this manner 
dissipation enters into our results. 

For the exactly solvable dissipative harmonic oscillator this programme could be carried 
through without any approximations. Within certain approximations the spin-boson model 
could also be mapped onto an unitarily equivalent spin-boson model by decoupling all except 
the nearly resonant bath modes. This procedure is depicted in Fig. 8. Large energy differences 
are decoupled for small flow parameters £ and smaller energy differences only later for larger 
flow parameters. This yields the fundamentally important separation of energy scales underlying 
e.g. renormalization theory. Notice the main difference from adiabatic renormalization where 
decoupling starts from the UV-cutoff and is never concerned with the low-lying states, also 
their energy difference to the discrete states might become comparable to the decoupled states 
at some stage of the renormalization procedure. Hence in order to obtain an effective spin-boson 
model with only nearly resonant couplings remaining, the flow equation technique is more suited. 
Using infinitesimal unitary transformations it is also apparent that the observables have to be 
transformed as well as is essential when states with an energy difference of order A r are beginning 
to be decoupled. 

As we have explained such a two-level model with only nearly resonant couplings becomes 
effectively equivalent to a dissipative harmonic oscillator as the higher states of the harmonic 
oscillator are not occupied. The exact solution of the dissipative harmonic oscillator could then 
be used as an asymptotic conserved quantity to calculate the correlation functions of the effective 
spin-boson model. This program could be carried out for an Ohmic bath with small coupling 
(a ;$ 1/4) and for super-Ohmic baths with arbitrary coupling. The main restriction was that 
the temperature should be much smaller than the renormalized tunneling frequency A r as only 
then the equivalence to the dissipative harmonic oscillator is possible. 

We have found agreement of our results with the NIBA on intermediate time scales set by A r 
where the NIBA can probably be trusted. The long-time behaviour of the spin-spin correlation 
function that is not accessible by the simple NIBA has an universal algebraic decay put forward 
by Sassetti and Weiss 18]. Their generalized Shiba-relation connecting intermediate and long- 
time scales was found to be fulfilled within numerical errors for super-Ohmic baths. Considering 
the fact that this generalized Shiba-relation is a highly nontrivial universal property far beyond 
a "simple" universal algebraic decay C(t) oc t _s_1 , this indicates that the approximations in this 
paper retain all the "relevant" physics. 

Compared to our approach the other methods mentioned in the Introduction have specific 
shortcomings in the low temperature regime. E.g. quantum Markov processes yield an entirely 
different long-time (low-energy) behaviour. There correlation functions typically show an expo- 
nential decay, even for the exactly solved harmonic oscillator. This can only be obtained from a 
Hamiltonian with a spectral function of the bath that is strongly modified with respect to its low 
energy properties as can easily be seen from the exact solution. However, a quantum Markov 
process can at least be derived from a Hamiltonian description of system plus bath since the 
complete positivity condition is satisfied pj|. For the NIBA already this point seems less clear. 
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At present it is unclear whether our method with some other approximations can also be used 
in the left -open parts of the parameter space of the spin-boson model. In particular it would be 
interesting to investigate temperatures T > A r . 
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Appendix 1 

In this appendix we briefly describe the adiabatic renormalization procedure applied to the model 
( |2.l| ) in the spirit of Wilson Q . In a renormalization procedure the Hamiltonian is divided into 
two parts H = Ho + Hi, where Ho describes the high-energy modes and Hi describes the 
low-energy modes. The Hamiltonian is then mapped to an effective Hamiltonian H^- 1 ' that 
operates on the Hilbert space of the low-energy modes and has, at least approximately, the 
same low-energy spectrum as H. If necessary this procedure can be iterated and one constructs 
a series of Hamiltonians H^. In the simplest perturbative renormalization scheme one sets 
= Po^Pq where Po is the projector onto the ground states of Hq. This is the first step in a 
degenerate perturbational treatment, higher orders can be calculated similarly. In the adiabatic 
renormalization scheme, the Hamiltonian Q2.1| ) is first transformed using a unitary transformation 

U = exp( V **A(b k - b{)) (A.l) 

where Iq = [uj,uj c ]- The unitary transformation does not affect A, but only Hs- The transformed 
Hamiltonian becomes 

= + £ (A.2) 
where 

H®=VZHsU -A* £ ^ (A.3) 

This Hamiltonian is now projected onto the ground states of Y2u> el ^k^fik- Thereby we obtain 
an effective Hamiltonian which has the same form as the original Hamiltonian, but with a new 
cutoff Q and with a new renormalized Hs- oj has to be large compared to typical excitation 
energies of Hs- The projection onto the ground states of Yju> ei u kb\bk corresponds to the first 
order of a degenerate perturbational treatment. Let P be the projector onto the ground states 
of Y2u k ei Wfc&fc&fc. The next term in the perturbational treatment is 

PTJlH s Uo(l - P){ u k bib k )-\l - P)W H S U P. (A.4) 

This term is of the order of a typical squared excitation energy of Hs divided by Co. It is small 
compared to the first term as long as u is large compared to typical excitation energies of Hs- 
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This estimate shows one of the advantages of the adiabatic renormalization scheme: Corrections 
can be estimated and are small as long as the effective cutoff is large. 

The procedure can now be iterated. We can introduce a new lu that is smaller than the old 
one but large compared to the typical excitation energies of the renormalized Hs- Alternatively, 
Co can be determined self-consistently to be large compared to typical excitation energies of the 
renormalized H$- The final result of this procedure is 

H cS = H S)CS + {AX k (b k + bl)+ LO k blh) (A.5) 

with 

H s , e s = ( e ' A ^ a< "k<^^^- b h Hse A ^^<^ k <^^( b ^ b l)^ B _ A 2 K ( A _ g ) 

U)<U) k <U) c 

{.) b denotes an average over the bath. This new Hamiltonian contains no high-energy scale. The 
initial Hamiltonian has to contain a counterterm of the form A 2 A^/u^. This is well-known 
in the theory of dissipative quantum systems. The main problem is that it is still difficult to 
calculate dynamical correlation functions using the renormalized Hamiltonian since it still has a 
complicated structure. 
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Figure captions 



Fig. 1. Testing the quality of the numerical routines for the exactly solvable dissipative harmonic 
oscillator. Here the spectral function is Drude-like J(u>) = y 2 oja/ '(y 2 +lu 2 ) with parameters 
a = 0.1,7 = 1.0 and Aq = 1.0. The function K(lo) defined in ( 3.46| ) determines all the 



equilibrium correlation functions and is plotted from the exact solution ( |3.50| ) and the 
numerical solution of the flow equations. Here A = 0.1 was chosen, but no dependence on A 
can be observed within numerical errors. 

Fig. 2. Results from the numerical solution of the flow equations for the spin-boson model with 
an Ohmic spectral function J(u>) = 2au)Q(u c — u). Parameters are a = 0.1, u c = 10.0 and 
Ao = 1.0. In Fig. 2a the flow equations are numerically integrated up to A = 0.2 and then 
the asymptotic conserved quantity is added to C(w,4). In Fig. 2b A = 0.1 is chosen and in 
Fig. 2c we have integrated even further until A = 0.05. The final result for the correlation 
function C(u>) is within small errors independent of these values A as can be seen in Fig. 2d. 

Fig. 3. Correlation functions C(u) for the spin-boson model with an Ohmic bath (u c = 10.0, Ao = 
1.0) and various values of a. The functions C(u) are plotted as functions of u/A r with 
the low-energy scale A r used as a fit parameter to identify the peaks. A r is set by ( |5.11 ). 



Fig. 4. Comparison of the flow equation result for the correlation function C{u) with the NIBA- 
curve. An Ohmic bath is used with parameters uj c = 10.0, a = 0.1 and Ao = 1.0. For 
intermediate time scales the NIBA-results agree with the flow equation results, however, 
the long-time behaviour of the NIBA is wrong as can be seen from the divergence of the 
response function S{u) for uj — ► 0. The data points for C(uS) show how the bath energies 
are discretized in the flow equation framework. 

Fig. 5. The normalized response function S(u) is plotted as a function of the rescaled energy uj/A r 
where A r is again used as a fit parameter to identify the peaks. An Ohmic bath is used 
with parameters a = 0.1, Ao = 1.0 and various values of uj c . Universality is confirmed with 
excellent accuracy. 

Fig. 6. Comparison of the flow equation result for the correlation function C(co) with the NIBA- 
result. An super-Ohmic bath with a spectral function like in Eq. ( |5.12| ) is used. Parameters 
are s = 2, K = 5, Ao = 1.0 and uj c = 10.0. The main difference to the NIBA is again the 
long-time behaviour as can be clearly seen in the log-log-plot. 



Fig. 7. Correlation functions C (w) for the spin-boson model with a super-Ohmic bath of type ( 5,12j ) 



Parameters are s = 2, Ao = 1.0, uj c = 10.0 and various values of K. A r is again used as a 
fit parameter to identify the peaks. It is set by the low-energy scale derived from adiabatic 
renormalization for s = 2: A r = cAo exp(— uo c /{2K)). c is a constant of order 1. 

Fig. 8. Sketch of the flow of the Hamiltonian and of a generic observable as a function of I in 
different regimes. For small values of I the system finds the low-energy scale. In this region 
the flow of the observables is negligible, and the flow equations are equivalent to adiabatic 
renormalization. In the other regions, the flow of the observable becomes important. The 
w-scale in the plots of J(u>,£)/J(lu,0) is logarithmic. 
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